function distance = euc_dist(x1,y1,x2,y2) %distance between 2 points (x1,y1) and (x2,y2)
    distance =sqrt((x2-x1)^2+(y2-y1)^2);
end    